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Abstract 


Conventional heavy fermion theories require existence of massive fermions. We 
show that heavy fermion phenomena can also be simply explained by existence 
of bosons with moderate mass but temperature dependent concentration below 
the formation temperature Ta, which in turn is close to room temperature. The 
bosons B ++ are proposed to be in chemical equilibrium with a system of holes h + : 
B ++ ^ h + +h + . This equilibrium is governed by a boson breaking function f(T), 
which determines the decreasing boson density and the increasing fermion density 
with increasing temperature. Since HF-compounds are hybridized from minimum 
two elements, we assume in addition existence of another fermion component h/ 
with temperature independent density. This spectator component is thought to 
be the main agent in binding the bosons in analogy with electronic or muonic 
molecules. Using a linear boson breaking function we can explain temperature 
dependence of the giant linear specific heat coefficient 7 (T) coming essentially 
from bosons. The maxima in resistivity, Hall coefficient and susceptibility are 
explained by boson localization effects due to the Wigner crystallization. The 
antiferromagnetic transitions in turn are explained by similar localization of the 
pairing fermion system when their density nh(TFL) becomes lower than nwc, the 
critical density of Wigner crystallization. The model applies irrespective whether 
a compound is superconducting or not. The same model explains the occurrence 
of low temperature antiferromagnetism also in high-T c superconductors. The 
double transition in UPta is proposed to be due to the transition of the pairing 
fermion liquid from spin polarized to unpolarized state. 


I. INTRODUCTION 

Despite of the great efforts to understand heavy fermion (HF) properties [1-3] in terms of 
conventional Fermi liquid theories [4-6], the Kondo model [7] or unconventional order param- 
eter models [8,9], the theories have not provided a clear overall picture of the most dramatic 
properties: the large value and the temperature dependence of specific heat coefficient 7 (T), 
the resistance maximum at low temperatures, the Hall coefficient and the susceptibility maxima 
near the same temperature. Likewise the superconducting states exhibit power law behaviors 
of several quantities such as the specific heat near T = 0 and an additional linear term in the 
specific heat. From the above reviews the following picture emerges: In the HF-compounds 
the /-electrons show localized behaviour at high temperatures and become delocalized at low 
temperatures. The point of division for the two behaviours is the Kondo temperature Tk • 

The purpose here is to show that the spectator fermion superfluid model (SFS), originally 
proposed for high-Tc compounds [10,11] and the superfluid states of heavy fermions [12,13], 
also provides a global picture which connects the superfluid state with the normal state and, 
at the same time, explains the properties of non-superconducting HF-compounds within the 
same framework with bosons decaying into pairing fermions, the holes. Since the boson density 
for T ^ T c and the pairing fermion density for T< T c both become small, they get local- 
ized due to Wigner crystallization [11,14]. In what follows we intend to show that the boson 
localization temperature Tbl < T# is in fact the coherence temperature corresponding to max- 
ima of Hall coefficient and susceptibility. Similarly the pairing fermion Wigner crystallization 
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temperature is Tfl = 7jv, the Neel temperature, since the fermion spins are antiferromagnet- 
ically ordered for T < Tfl • Both temperatures Tql and Tfl are logical consequences of the 
SFS- model, not an extra input. In the case of high-T c compounds of the type RBa 2 Cu 307 _ x , 
3D-antiferromagnetism is known to exist together with superconductivity [15,16], with the Neel 
temperatures in the range 0.5 K< Ts <2.2 K. For HF-compound UPt 3 the low temperature 
peak at Tpl « 20 mK observed by Schuberth ei ah [17] is here proposed to be also an example 
of Wigner crystallization of the pairing holes. Immediately above Tfl the holes are proposed 
to be polarized up to the lower temperature T™ of the double peak of UPt 3 as a consequence 
of the ordering of the pairing hole liquid, rather than having unconventional order parameter 
for the superfluid. The sequence of temperatures Tfl = T Xr, T~ and T c + is in agreement with 
MC-calculations by Ceperley et al [18,19]. 

The plan of the paper is as follows. In chapter II we propose that the giant linear specific 
heat coefficient 7 (T) is a consequence of the bosons as charge carriers in these compounds. 
In chapter III we show how the maxima in Hall coefficient and susceptibility follow from the 
boson localization. Knowing the boson localization temperature Tbl > T c , one can estimate 
the pairing fermion localization temperature Tfl = TV, which is discussed in chapter IV. 
The Wigner crystal has two transverse electronic sound modes, and in the fermion localization 
case also spin waves with a gap. The longitudinal plasmon of fermion Wigner crystal gives a 
spectrum with temperature dependent gap below Tn , similar to the spin wave spectrum. 


II. THE NORMAL STATE SPECIFIC HEAT 

We follow closely the discussion of SFS carried recently [14] through in details for the high-T c 
compound YBa 2 Cu 3 _ x M x 07 _$ and proceed to show that the same model can explain also many 
of the heavy fermion features mentioned above. We adopt here the same simple principle proved 
successful in the case of 123-compounds: The temperature dependence of an experimental 
quantity is simply determined by the concentration dependence of that same quantity, since 
carrier densities are temperature dependent according to the charge conservation in the chemical 
equilibrium 


n B (T) = n 0 f(T) 

n*(T) = 2no[l-/(T)] (1) 

n t (T) = constant, 

with /( 0) = 1 and /(7 b) = 0. These equations are quite general and form a cornerstone for 
SFS. In particular we don’t need to assume statistical mechanics of non-interacting particles 
to their validity. The holes with temperature dependent density n/j(T) are treated in the 
superfluid and normal states like conventional charge carriers. Their density of states can 
in principle be calculated from the two-fluid formalism [14]. The fact that their density of 
states is temperature dependent can be looked upon as coming from the existence of bosons. 
With the concentration rule we may proceed to connect various experiments without knowing 
the exact band structure and the binding mechanism of the bosons [11,14], In the normal 
states we use a simple linear function f(T) = 1 — T/T B . Since T c ’s are zero or very small 
in comparison with the boson formation temperature T B we can take f(T c ) « /( 0) = 1 , as 
a first approximation. These features come out of an explicit fit of the model parameters to 
experiments in the case of UPt 3 and URu 2 Si 2 , to be discussed later in more detail. We assume 
the existence of bosons also in the non-superconducting cases. The picture proposed here is 
very close to ionization of atoms in gas. There the fraction of non-ionized atoms is given [20] 
by f(T) = 1 — [1 -f (P/Po)(7b/T) 5 ^ 2 e 2 ^ fcBT ]“ 1 ^ 2 , which is linear in a range of temperatures. 

Following the concentration rule, the internal energy is taken to be a sum of a boson contri- 
bution and the pairing and the spectator fermion contributions 

u(T) = |fc fl n 0 T/(T) + i Th T 2 [l - f(T)) + | T ,T 2 , (2) 

where 7 * and 7 * are linear specific heat coefficients at T = 0 for the two fermion components. We 
will use the expression for noninteracting bosons, and therefore Eq. (2) represents the simplest 
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approximation one can think of. For / = 0 this reduces to the conventional two fermion result 
and for / = 1 we obtain the free boson and the spectator contributions. In the normal state this 
expression is supposed to be valid only down to finite temperatures near T c where the superfluid 
fluctuations become important. In the superconducting state one may use the expression (2), 
except that the boson part is replaced by the acoustic plasmon contribution to be discussed 
later. In the BCS-like situation ( Tb = T C1 no plasmon sound) we have f(T) = n < (7 1 ), where the 
superfluid fraction behaves exponentially, n, « 1 — with t = T/T c , Eq.(2) reproduces 

the typical BCS result with a step in the specific heat at T c . The expression (2) is therefore 
reasonable in the three special cases / = 1, / = 0 and / = n,. We anticipate that the calculation 
of boson breaking function f(T) is in principle possible once we know the details of the boson 
binding, the band structure and the chemical lattice. So far we are content to assume that the 
spectators play an essential role in the binding [14]. In view of the fact that the compounds we 
are working with are rather complicated, we therefore should use a formalism, simple enough, 
to enable one to relate various experiments with a reasonable accuracy as was done in Ref. [14] 
for the 123-case. In what follows we intend to show that also the heavy fermions can be treated 
in this spirit and the results turn out to be equally good. 

Eq.(2) gives for the normal state the specific heat coefficient 


7 (T) = 


c(T) _ 3 feanp 
T “ 2 T ' 




(3) 


using the simple linear function for /(T). At low temperatures the boson contribution domi- 
nates, since Tb ~ 300 - 500 K, and at low temperatures ( T < 20 K) we obtain 

7 (T) = 7. + (4) 

This also shows, how the spectators give the linear term in the specific heat in the supercon- 
ducting state at T = 0. At the temperature Tb the linear term reads 

7CO=7.+p». (5) 


Above Tb we obtain the normal metallic behaviour 7 (T) = 7, + 7*. Immediately below Tb 
we obtain linear behaviour and a step at Tb- We anticipate that mj >> mj, therefore the 
contribution from the pairing fermions 7 a should dominate at room temperature. On the other 
hand, already at much lower temperature the boson contribution gets small as compared with 
fermion parts, therefore one may be able to separate the contributions from 7* and 7, , since the 
latter one can be obtained at least in the superconducting situation from the measurements near 
T = 0. Near the localization temperatures Tbl and Tfl these expressions have to be modified 
in a way to be discussed later. In particular the boson term cannot continue down to zero 
temperature but is intruded by the superconducting or the antiferromagnetic transition. Since 
these are ordered states, the boson contribution is replaced below the ordering temperature by 
the plasmon sound or spin wave contributions, which will give finite entropy. 

We have compared the electronic specific heat coefficient y(T) from Eq. (4) with experiments 
in Fig. 1 a) for UBei3 and CeCu 2 Si2 and in Fig. 1 b) for CeAU and CeCu 6 . These figures are 
based on figures (6) and (16) of Stewart [2]. It is seen that this simple theory gives a satisfactory 
explanation for the observed behaviour of 7(T) in these four cases. The coefficients A come 
out to be nearly the same A « 1400 mJ/Kmol corresponding to boson densities of the order 
10 21 cm~ 3 . One obtains this order of magnitude also from the the analysis of the Hall-data, 
as we shall see. The model therefore proposes that the large size of 7(T) is due to the bosons 
rather than fermions. This means also that the fermions need not be excessively heavy, which is 
desirable, since band structure calculations do not produce excessively narrow bands required by 
the conventional heavy fermion theories (see Ref. [4], p. 121). In fact the bandwidths obtained 
from band calculation are two orders of magnitude larger. 

The compounds in Fig. 1 a) become superconducting at very low temperatures, but those in 
Fig. 1 b) do not. This simple theory should apply to both types. As we shall see the present 
model predicts a glassy antiferromagnetic transition at very low temperatures so that 7(T) 
is cut off. For CeAb this is clearly visible (see Ref. [2], Fig. 15) but for CeCu6 less so. At 
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present we are not able to determine the fermion effective masses very accurately. There is 
some uncertainty associated with the possible non-linearities in /(T) which no doubt can be 
estimated easily, but goes beyond the scope of this paper. Clearly Eqs. (2)-(5) cannot be valid 
down to zero temperature, as mentioned above. For superconducting compounds 7(T) near T c 
is modified by superconducting fluctuations and below T c the plasmonic contribution replaces 
the free boson part. For non-superconducting compounds the superconducting fluctuations may 
exist near T = 0. The ordered transitions at T e or Tfl = Tjy seem to offer a natural cut off for 
most HF-compounds. 
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FIG. 1. The electronic specific heat coefficient y(T) for a) UBeia and CeCu 2 Si2 and b) CeAhj and 
CeCu 6 . Solid lines presents our theoretical result from Eq. (4) with parameters: a) A = 1300m J/molK 
and 7, 2* 5mJ/molK 2 and b) A = 1500m J/molK and y g 5* 140mJ/molK 2 . Experimental points are 
from Figs. 6 and 16 in Ref. [2]. 


III. BOSON LOCALIZATION 

A logical consequence of the SFS model is the Wigner crystallization (WC) of bosons, because 
at a certain temperature Tbl, which lies in between T c and Tb 9 their density diminishes below 
the critical value corresponding to r s = 170 [19]. The localization temperature Tbl can be 
estimated from the condition for the critical density, 

nwc = n 0 f {Tbl) « n 0 (l - T BL /T B ) « 0.16 • 10 21 (^) 3 cm- 3 , (6) 

where m is the electron mass and the boson mass m# = 2m^. This order of magnitude 
estimate corresponds to a jellium background and may be inaccurate for backgrounds prevalent 
in the heavy fermions with spectators and the chemical lattice. Most likely, however, such 
a transition does take place, and the main consequence is that the bosons become gradually 
delocalized below Tbl • The measured carrier densities in HF-compounds are of the order of 
magnitude 10 21 cm“ 3 [6], close to those required by Eq. (6). Above Tbl the bosons cease to be 
active charge carriers. 

A fully consistent thermodynamical description of the localization is not attempted here for 
the following reasons: Above Tbl the bosons continue to decay into holes with approximately 
linear temperature dependence, as deduced from Hall-measurements. Because of this decay, 
the Wigner crystal is never complete. Since the boson localization phenomena occur at mod- 
erately low temperatures, the Wigner crystal is never in thermal equilibrium, and one expects 
a glassy behaviour. The situation is further complicated by the lattice which can cause par- 
tial localization of the bosons even below Tbl • This is presently the main uncertainty which 
clearly affects our calculation of the coefficient A in Eq. (4) with density no deduced from the 
Hall-measurements. 

The Monte-Carlo-calculations by Ceperley et al. [18,19] indicate, that the energy densities of 
Coulomb liquid and crystal states are not very different, therefore the associated specific heat 
jump may be small at elevated temperatures. For these reasons the liquid- crystal transition is 
very rounded and shows up most clearly in the Hall density which well above Tbl is determined 
solely by the holes. This explains the near linearity of n# (T 1 ) observed in virtually all compounds 
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above the coherence temperature, which we associate here with Tbl . We actually do not 
exactly know what /(T) to use between T c and T B L t therefore we have used the following 
purely phenomenological description: Divide the density of bosons below Tbl into localized 
and delocalized portions with a phenomenological function £(T) such that the density of mobile 
bosons is 


n B DL(T) = n 0 /(T)€(T) « n 0 (l - T/T B ) £(T), (7) 

where ^ = a < 1 for T « Tbl and £ = 0 for T » Tbl • Here a < 1 means that a fraction 
(1 — or) of the bosons remain localized for T « T B l . Localized bosons continue to produce 
holes with linear rate /(T). This was previously [14] used successfully to explain the Hall effect 
in the high-T c case. The existence of bosons in the normal state was also used to explain the 
"reservoir” effect for 123-compounds. 

The function £(T) can be obtained from the experimental Hall coefficient R a b oc n^(T)~ x . 
The effective Hall density (in the a6-plane) reads 

_ f 2n 0 (l - f(T)) =* 2n„£, for T > T Bt . . 

\ 2n„(l - /(T)) + 2n 0 /(T)£(T) 2 2 n„£ + 2n 0 (l - T/T B )((T), for T <T BL . U 

These equations together with the Hall-measurements can be used to fix T B as well. The 
localization explains the minimum in n#, a feature which seems to be as common for heavy 
fermions as it is for high-T c superconductors [14]. We have shown for 123 compounds [14] that 
f(T) is a universal function. 

The Hall density in the normal state together with the experimental data of Lapierre et 
al. [21] and Schoenes et al. [22] is shown in Figs. 2 a) and 2 b). The theoretical curve in Fig. 2 
a) for UPt 3 was calculated using Eqs. (8) with the phenomenological boson localization factor 
of the type 


t(T) = exp [3.5 ( ^ok )] • ' 6 K < T < 50 K ‘ 
Similar function was used also for URujSb in Fig. 2 b). 
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FIG. 2. a) The Hall density in the normal state for UPt 3 from Eq. (8) together with the experi- 
mental data of Lapierre ef al. [21]. Theoretical result is obtained by using parameters: Tb = 500 K, 
fc = .88, n 0 = 4.95 • 10 21 cm” 3 and a = 1. b) The Hall density in the normal state for URu^Sia 
from Eq. (8) together with the experimental data of Schoenes et al. [22]. Theoretical result is obtained 
by using parameters: Tb = 350 K, f c = .98 and no = 5.7 ■ 10 21 cm~ 3 and boson localization factor 


£(T) = .36 • exp 



for 18 K < T < 110 K. 


We overlook the further complication, that the spectators could bring in a more complicated 
two-band Hall behaviour. The fact that, Hall coefficient with H ± c is nearly constant whereas 
the one with H || c changes with temperature above T B l [23], gives evidence that band mixing 
is small, therefore this point is not fully understood at present. The outcome of Eq. (8) for 
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two compounds, UPt 3 in Fig. 2 a) and URu 2 Si 2 in Fig. 2 b) gives us their 7js-values and 
the localization temperatures. The reason for selecting these two HF-compounds is that they 
represent the two main cases with respect to fermion localization: In UPt 3 the antiferromagnetic 
transition turns out to occur in the superconducting state and for URu 2 Si 2 in the normal state. 
In both cases n/f-curves show slight upwards bending for T >Tbl and the Tb- values are of the 
order 300-500 K. The localization temperatures are Tbl & 30 K for UPt 3 and approximately 
50 K for UR^Sh- 

Also the experimental specific heats show structure associated with Tbl and Tb. The ap- 
proximate separation of the electronic part AC by Renker et al [24] shows a broad maximum 
for UPt 3 near 30-40 K and a peak near 300-500 K. We expect that the pairing fermion band 
and the spectator band rearrange themselves near Tb, which could also show up in the elec- 
tronic specific heat. The broad maximum at 30-40 K can be interpreted in terms of Bose glass 
transition associated with boson localization. A qualitatively correct behaviour is obtained by 
taking a linear superposition of the Bose liquid and the crystal internal energies \ns{T)kBT 
and 3 ub(T)A:bT for the boson liquid and solid, respectively, hence 


A C = 


f "o k B 4, [T/(T)( 2 - ((D)] + y.T + §7h [1 - /(Dl T , for T < T B 
(is + lh)T , for T > Tb 


( 10 ) 


The total electronic specific heat from Eqs. (10) is compared with the data of Renker et al. [24] 
in Fig. 3 for UPt 3 . The agreement is rather striking. In particular the boson formation temper- 
ature Tb shows up clearly in both compounds as a peak with quadratic temperature dependence 
near Tb. Just above Tb we predict a jump by amount J IhTb in the specific heat. Unfortunately 
the data do not go high enough in temperature to show the jump. In the case of UR^Sh one 
obtains from Renker et al. [24] similar picture except for the peak near the antiferromagnetic 
transition at Tat = 18 K, which will be discussed later. In the latter case the jump at Tb shows 
up. 



Temperature (K) 

FIG. 3. The electronic specific heat for UPt 3 after Renker et al. [24]. Solid line is our result 
from Eq. (10 with parameters: y, = 0 and y h = 16mJ/molK 2 . Dash-dotted line is the result with 
y» = 4mJ/molK 2 and y h = 16mJ/molK 2 . 

The boson localization temperature Tbl limits the attainable transition temperature T c , 
since Tbl > T c . By Eq. (6) one should make the effective mass m* small to increase Tbl an d 
thereby the T c . At the same time the background can also have a decisive effect. In Ref. [14] 
we found effective mass ~ m e and the estimate for UPt 3 is 5-20 electron masses. 


IV. FERMION LOCALIZATION AND THE ANTIFERROMAGNETISM 

Since the Wigner crystallization takes place independent of the statistics, also the pairing 
holes must become localized at some temperature Tfl < Tbl , since their density vanishes at 
T = 0. Assuming that Tbl is known, we can actually calculate Tfl from Eqs. (1) since the 
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bosons have double charge and mass. The crystallization densities of holes and bosons are 
therefore related by 


2n 0 [l-/(7> L )] = 


npf(T B L ) 0 
64 W 


( 11 ) 


The factor Q takes into account the fact that the background here is not jellium, which would 
correspond to Q = 1 and r, = 170. For the superconducting compounds we have to make 
distinction between the two cases Tfl <T c and Tfl > T c . It is well known that for fermions the 
spins are in antiferromagnetic arrangement in the Wigner crystal phase for jellium background. 
In this connection one should remember that for somewhat higher densities, hence here above 
Tfli the holes may form a spin polarized liquid. The direction of spins may be strongly 
influenced by non-jellium background effects, but WC-localization offers a natural explanation 
for occurrence of antiferromagnetism in the HF-compounds. If the antiferromagnetic transition 
takes pk e above T c we obtain from Eq. (11), by using the linear approximation for f(T)> 


Tfl = — (1 - — -)Q « —Q. ( 12 ) 

FL 128' T b )W 128 V 1 ’ 

Depending upon the values of Tb and Q this gives Tfl > 1 K. The main source of inaccuracy 
in Eq. (11) is the factor Q and the sensitivity of the left hand side to possible non-linearities in 
/(T), which show up in the experiments for the Hall densities illustrated in Figs. 2 a) and 2 b). 

The fermion localization at Tw = 18 K explains also the sudden drop in for UR^Sh 
apparent in Fig. 2 b). The amount of the drop is larger than we would expect. This we 
believe is coming from spectators and for a quantitative calculation of the drop one should 
use a two-band formula. The main point we want to make is that the strange behaviour of 
URu2Si2 Hall coefficient is understood with the boson and fermion localizations. This includes 
also the different behaviour for UPt3. Conversely the Hall density is like a map where the 
localization effects can be read out. 

To carry out a calculation of Tfl one must first determine f(T) from experimental data on 
Hall effect, specific heat etc.. As an example we show how to calculate Tfl for UPt3, where 
the antiferromagnetic transition supposedly occurs in the superfluid state. Another example 
we have been able to carry through with reasonable accuracy is the superconductor URu2Si2, 
where Tn — Tfl ~ 18 K in the normal state. We start the discussion with UPt3. Using the 
Hall data, the specific heat data and the exponent relations derived by Kallio et al [12] we 
have determined the boson breaking function f(T) and the superfluid fraction n,(T) using the 
functional form ( t = T/T c ) 


/(t) = n,(0 + /(T e )< 4 - 5 ^ 3 

n,(t) = 1 — at 7 — (1 — (13) 

The philosophy behind these terms is of course that the functions f(t) and «s(0 have direct 
connection with experiments and the exponent 7 is limited by the two-fluid model to the range 
0 < 7 < 1.5. For UPt 3 we use the specific heat data of Midgley et al. [25] and Hall density data 
of Lapierre et al. [21] to obtain for the parameters the values 7 = 1.42, f(T) = .88, a = .91, 
and k = 3.5. In the normal state we therefore obtain 

f(T) = .88(1 - T/T b ) ,T b 2 500 K. (14) 

In calculating the specific heat we used Eq. (2) with the boson contribution replaced by the 
plasmon term [14] 


^ (k B T? 

Upl ~ 30[/i«i(T)] 3 ’ (15) 

where ui(T) is the plasmon sound velocity (see Eqs. (16) and (17)). The data for UPt 3 are 
extremely well fitted with these functions. As an example we show the specific heat in the 
superfluid state in Fig. 4 a). Because of the low value of T C} f(T c ) is nearly unity justifying our 
approximation f(T) = 1 -T/T B in Eq. (3). The value of exponent 7 = 1.42 in Eq. (13) is close 
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to the one obtained for 123 in Ref. [14] (7 = 1.3) as well as the gap parameter k = 3.5 instead 
of 3. The small coefficient .09 in front of the gap term makes n 9 (t) to have nearly power-law 
behaviour. The most dramatic difference between the HF- and high-T c compounds is the pair 
breaking at T e which gives f c « 1 for former and f e « .6 for the latter. For URu 2 Si 2 it comes out 
close f e ^ .98. The physical explanation for this is the low value of T c for HF-compounds. The 
heavy fermions are therefore further away from the BCS than the high-T c superconductors, in 
view of the fact that within the present framework f c = 0 for BCS. Since the boson densities are 
not very different for the high-T c and HF-compounds, simple Bose-Einstein model for heavy 
fermions would give, with moderate boson mass ms < 10m e , right order of magnitude for 
T c [14,26]. Without boson breaking this model would, however, have difficulties in explaining 
the localization phenomena discussed here. 



FIG. 4. a) The experimental specific heat for UPt 3 from Ref. [25] together with our theoretical 
result. For f(t) we have used parameters: f c = .88, 7 = 1.42, a = .91 and k = 3.5. The used spectator 
contribution 7, = 130mJ/molK 2 is much larger than the one obtained from Fig. 3. b) The experimental 
specific heat for URu 2 Si 2 from Ref. [30] together with our theoretical result. For f(t) we have used 
parameters: f c = .98, 7 = 2.16, a = .987 and k = 3.6. The spectator contribution comes out to be 
7 j = 23mJ/molK 2 . 


Knowing the boson breaking function f(T ) for UPt 3 we can now calculate the Wigner crys- 
tallization temperature Tfl for the holes to be about 20 mK, close to the experimental value 
18 mK found for the extra peak by Schuberth et al. [17]. Our result is obtained from Eq. (11) 
with Q=l. Since in the Wigner crystal for the holes the spins are in the antiferromagnetic 
arrangement, we propose that Tfl = TV, the Neel temperature. 

Although the properties of the antiferromagnetic transition may be affected by the chemical 
lattice, we may explain some properties, like the specific heat peak, by the Wigner crystal 
model below TV. Since the transition takes place at low temperatures, we may use the sound 
mode approximation for the WC specific heat. The excitation modes of WC have been studied 
theoretically by Carr [27]. According to his calculations the Wigner crystal has two transversal 

sound modes, with the density dependent sound velocities oc rf 1 ^ 2 . In the present case this 
translates into the temperature dependence [14] 

(16) 

In the case of pairing fermions forming a liquid, the longitudinal sound velocity ui(T) was 
determined by the relation 

(17) 

Despite of the pairing fermion localization the two-fluid philosophy remains the same: The 
normal liquid density in the superfluid system consisting of bosons and holes, is obtained from 
the bosonic (sound modes, spin waves) and fermionic excitations (boson breaking). The latter 
ones are determined by [1 — f(T )] which should be determined by the chemical equilibrium, but 
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we choose to follow the philosophy of Ref. [14] and try to determine f(T) from experiments, 
if possible. One uncertainty remains: Are the sound mode or the spin waves the dominating 
feature in the bosonic excitations, which determine the specific heat in the superfluid state? 
The data on UPt 3 is insufficient to give an answer, because Tfl is so low. If the transition 
takes place in the normal state, we believe that immediately below Tfl the boson part in the 
specific heat will be dominated by the spin waves as compared with plasmonic contributions. 
This will be apparent in URu 2 Si 2 , to be discussed later. Much below Tfl the spin wave 
dominated specific heat would go to zero exponentially when T — ► 0 since there is the gap, 
which experimentally turns out to be rather large. We expect this to be true also for non- 
superconducting compounds: The giant bosonic specific heat coefficient is cut off at Tfl = TV , 
because in the antiferromagnetic state the bosonic coefficient has to be replaced by the plasmonic 
contribution. We of course do realize, that some of the HF-compounds, so far, show no order 
whatsoever, below the lowest temperatures measured. One such example is CeAl 3 , with 7(0) « 
18.5 mJK -2 cm” 3 . However, below 0.32 K for this compound the specific heat coefficient 7 (T) 
starts reducing, which can be interpreted as coming from a glassy antiferromagnetic transition 
(See Ref. [2], Fig. 15). One should also beware of that some of the antiferromagnetic transitions 
can have a different origin. 

The antiferromagnetic Wigner crystal made of the holes is also never in equilibrium because 
their number density is changing with temperature and one has vacancies and/or the full WC- 
order can never be extended in practice. For this reason one should expect also here spin glass 
behaviour. In the high-T c case the situation should be similar: The antiferromagnetic transition 
of holes should take place at about 1-2 K, if one uses Eq. (12) and the f(T) determined in 
Ref. [14]. This is born out in the experiments on Bi 2 Sr 2 CaCu 2 0s by Caspary et ai [28], 
who separated out the electronic specific heat in this temperature range and found broad 
maxima of 7(T) at different magnetic fields. These maxima were interpreted as spin glass 
behaviour by Caspary et a/., for different reasons though. Several other high-T c compounds 
show antiferromagnetic transitions at low temperatures below T c . 

A further feature easily understood within the present model is that the antiferromagnetic 
specific heat peak should exist even in magnetic fields higher than /7 c2> because the bosons 
still exist, and hence their decay into holes occurs independent of superconductivity. This was 
indeed observed in UPt 3 by Schuberth et ai [29]. 



FIG. 5. The experimental specific heat for UPt 3 from Ref. [25] together with our theoretical result, 
with double peak at T c . For spin polarized liquid we have used parameters: T , T = -47 K, which 
corresponds T c “ = .904 • 7^ and k~ = 3.5. For unpolarized liquid: 7? = .52 K, k + = 4.3. Exponent 
7 = 1.41 and coefficient a = .91 are the same for both liquids. The gray line shows our result for 
unpolarized liquid below T c “. 


In the Monte Carlo calculations of Ceperley et ai [18,19] it was found that for r t < 100 — 170 
(hence T > Tfl in our case) one obtains a spin polarized liquid and below r, = 60 — 75 
unpolarized liquid. Since the energy difference between the two liquids is very small and in 
UPt 3 we have chemical lattice different from jellium, the possibility exists that the double 
peak near T e « .5 K is explained by the phase transition between the two types of pairing 
fermion liquids. In fact, if one makes use of f(T) deduced from experiments and assigns for 
T F l = .02 K the value r t = 170, we obtain at T c the value r, = 60, which is close to the 
value 75 obtained by Ceperley et ai for the crossing of the polarized and unpolarized liquid 
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energies. This interpretation of the double peak would be in agreement with the domain 
structure found by Midgley et ai [25], since the energies for the two structures are close. The 
associated domain structure could also be a reason for the metamagnetism observed in several 
compounds. Since the TVs are low it would also occur in the normal states. We stress that the 
present interpretation preserves 1 5o-order parameter for the superfluid. Actually the rvlimits 
for Wigner crystal of fermions are numerically uncertain even for jellium background due to 
the fact that the GMC-method is strictly speaking ” exact” only for Bose systems, as discussed 
by Ceperley et ai Since we have a different background, the right hand sides of Eqs. (11) and 

(12) can have sizable Q-factors, which can only be determined experimentally, owing to the 
complicated background and the glassy nature of the Wigner crystal transitions. We believe, 
however, that the order of the two localization temperatures remains always the same, i.e., 
Tbl > T n . 


V. CONCLUSIONS 

If we identify the pairing holes belonging to /-band in the SFS case the following picture 
emerges: The /-electrons combine into bosons at high temperatures T < Tb which remain 
localized down to temperature Tbl* In the temperature range Tbl < T < Tb the chemical 
equilibrium requires existence of mobile /-holes with density n* ~ 1 — /(T) « T/Tb , where 
/(T) is a universal function. At very low temperatures the pairing fermions suffer a glassy 
antiferromagnetic transition at Tjy. This explains the existence of superheavy quasiparticles in 
dHvA-experiments, which then, near Tn, measure the distribution of masses near a localization 
transition. 
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